Excitation spectrum of two correlated electrons in a lateral quantum dot 

with negligible Zeeman splitting 
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The excitation spectrum of a two-electron quantum dot is investigated by tunneling spectroscopy 
in conjuction with theoretical calculations. The dot made from a material with negligible Zeeman 
splitting has a moderate spatial anisotropy leading to a splitting of the two lowest triplet states 
at zero magnetic field. In addition to the well-known triplet excitation at zero magnetic field, two 
additional excited states are found at finite magnetic field. The lower one is identified as the second 
excited singlet state on the basis of an avoided crossing with the first excited singlet state at finite 
fields. The measured spectra are in remarkable agreement with exact diagonalization calculations. 
The results prove the significance of electron correlations and suggest the formation of a state with 
Wigner- molecular properties at low magnetic fields. 



Two-electron (2e) systems (e.g., the He atom I] and 
the H2 molecule |2J) have played an important role in the 
development of the theory of many-body effects in quan- 
tum mechanics and the understanding of the chemical 
bond. Recently, two-dimensional (2D) man-made quan- 
tum dots (QDs), fabricated in semiconductor materials 
have attracted attention as a laboratory for investiga- 
tions of few-body systems with highly controlled param- 
eters 0,SHI1,HIEI1,[I3] Such aS the electron number, 
confinement strength, interelectron repulsion, and the in- 
fluence of an applied magnetic field (B). Furthermore, 2e 
QDs may be used for implementing logic gates in quan- 
tum computing 

Earlier experiments on vertical QDs led to the observa- 
tion of a singlet-triplet (ST) transition in the 2e ground 
state as a function of the magnetic field 0, 0, an d f° r 
N > 2 to an energy level shell structure reminiscent of 
the periodic table of natural atoms [|| . The spin-triplet 
excited state of the 2e system was investigated in Ref.|7]in 
finite bias tunneling experiments. Laterally defined few- 
electron dots [8| added investigations of the magnetic- 
field dependent ST gap for two electrons 0, ^} . 

At low magnetic fields, the importance of electron 
correlations in 2D QDs (under appropriate conditions, 
see below), with different behavior compared to natu- 
ral atoms, has been anticipated theoretically [I2I fisl HiL 
Uti Firl llij . Clear experimental evidence for such 
strong correlation effects is still lacking. Insights can be 
gained through investigations of the excitation spectrum 
of 2e QDs in a perpendicular magnetic field. Here we 
report on joined experimental and theoretical endeavors 
that establish the correlated nature of the two electrons 
in our lateral QD, purposely fabricated to exhibit a neg- 
ligible Zeeman splitting. We find that the explanation 
of the measured spectra requires a high level treatment 
of electron correlations in the dot, and cannot be ac- 
counted for by perturbative schemes starting from the 
independent-particle picture. The exact-diagonalization 
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FIG. 1: (Color online) Differential conductance d//dVbias at 
zero magnetic field. Electron numbers N are indicated in the 
diamonds. The arrows mark the transitions between the one- 
electron ground state and the 2e spin-singlet (So) and spin- 
triplet state (Ti). Inset: scanning electron micrograph of the 
sample. The quantum dot is in marked with a dotted circle 
and the gate PG is the plunger gate. The vertical dashed line 
marks the bias voltage for which the data in Fig. [5] is taken. 



(EXD) calculations are in remarkable agreement with in- 
tricate measured features of the spectra, indicating the 
formation of an H2-type Wigner molecule (WM) at low 
magnetic field, lifting of degeneracies, and the appear- 
ance of avoided crossings between excited singlet states 
associated with a shape anisotropy of our dot - as well 
as excitations resulting from relative and center-of-mass 
(CM) motion. 

The sample is based on a 55 nm wide Ga[Al]As 
parabolic quantum well |l9| with the center located 
75 nm below the surface. The parabola ranges from 6% 
Al-concentration in the center to 23% at the edges. It 
was designed for having an average ^-factor for electrons 
\g\~0 20]. The effective mass of the electrons is higher 
than that of pure GaAs due to the larger band gap. We 
use the estimated value m* = 0.07toq, leading to an effec- 
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tive Rydberg energy — 6 meV. The two-dimensional 
electron gas (2DEG) has a density n B = 4.5 x 10 15 raT 2 
and a mobility n — 2.5 m 2 /Vs at a temperature of 4.2 K. 
A highly Si-doped layer 1.76 /im below the surface act- 
ing as a backgate was not used during the experiments 
described below. The quantum dot (geometric diame- 
ter 220 nm) with integrated charge readout depicted in 
Fig. ^ (inset) was defined using electron beam lithogra- 
phy patterned split-gates. All measurements were per- 
formed at an electronic temperature below 300 mK as 
estimated from the width of conductance resonances in 
the Coulomb-blockade regime. 

Figure ^ shows the differential conductance of the QD 
in the plane of bias and plunger-gate voltages. It was 
measured by applying the DC voltage V4ias/2 to the 
source and — Vbi as /2 to the drain contact. The recorded 
current / was numerically differentiated in order to ob- 
tain the differential conductance dl /dVbias- The resulting 
diamond pattern of suppressed conductance corresponds 
to fixed electron numbers N in the dot as indicated in 
the figure. 

The values for the electron numbers were found by 
depleting the dot down to zero electrons with more nega- 
tive gate voltages (not shown) and by using the quantum 
point contact near the quantum dot as a charge detector. 
For these modified gate voltage settings we observed the 
first single-particle excited state of the N = 1 quantum 
dot at an energy Ai = 5 meV above the ground state and 
found an addition energy A/i2 = 6.9 meV (both Ai and 
A/i2 are outside the range shown in Fig. 1). These num- 
bers show that in this dot, confinement effects (Ai) and 
interaction effects (A//2) are comparable in magnitude. 

Returning to Fig. ^ we determine from the extent of 
the N — 2 diamond in bias direction the addition en- 
ergy A/13 = 5.1 meV for adding the third electron to the 
dot. The strong lines labeled So in the figure correspond 
to tunneling transitions between the one-electron ground 
state and the 2e spin-singlet ground state. 

Strong lines parallel to the diamond boundaries outside 
the diamonds correspond to excited states of the level 
spectrum. In Fig.^ we see the transition T + between the 
one-electron ground state and the first excited 2e triplet 
state. The energy separation between the singlet and this 
triplet state is J = 1.7 meV. For negative bias voltage we 
see an additional excited state T_ which we attribute to 
the second triplet state split from T + by 1.5 meV due to 
the deviation of the dot potential from circular shape (see 
also Fig.0]below). Negative differential conductance can 
arise due to asymmetric coupling of states to source and 
drain. Inside the diamonds, higher order (cotunncling) 
processes can be seen for N > 2 21] . 

In order to investigate the evolution of the excited state 
spectrum with magnetic field, we have measured the cur- 
rent at fixed Vbias = 2.5 mV. In Fig.|3we plot the deriva- 
tive dl / dVpg. Resonances in dI/dV pg correspond to res- 
onances in the differential conductance dl/dVhi^s- The 
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FIG. 2: (Color online) Differentiated current dI/dV pe at 
Vbias = 2.5 mV. Red regions correspond to negative dI/dV ps , 
black to positive values. Electron numbers -/V are indicated. 
Transitions between the one-electron ground state and the 2e 
spin-singlet ground state (So), spin-triplet excited state (T+) 
and spin-singlet excited state (S2) are labeled. 

resonances indicated by symbols in the plot are labeled 
consistently with Fig. ^ R can be seen that the ST sep- 
aration J(B) — Et + — Es between the So and the Te- 
states decreases with magnetic field and becomes zero at 
Bst «4T. Beyond this field, the energy of the So state 
continues to increase and exhibits an avoided crossing 
with another excited state S2 which appears above 2.7 T. 
The avoided crossing suggests that S2 is an excited spin- 
singlet state and that the deviation of the dot potential 
from axial symmetry causes level repulsion. A third ex- 
cited state (T +j cm) appears at even higher plunger gate 
voltage. This state is shown below [see Fig. 4] to be a 
triplet excited state associated with the CM motion of 
the 2e molecule. 

A qualitative indication of the importance of correla- 
tion effects can be obtained via the Wigner-parameter 
i?w = (e / kIo)/^JjJq, i- e -j the ratio of the Coulomb in- 
teraction for two electrons at distance Iq = ^Jh/m*uio 
and the single-particle level spacing in a circular har- 
monic confinement, fuoo Correlation effects are of 
increasing importance when the interaction energy be- 
comes dominant over the single-particle spacing, i.e., 
when i?w > 1 . Using a conservative estimate for the con- 
finement energy of 5 meV (as quoted above for the dot 
with zero or one electrons), we find i?w ~ 1-55 indicating 
that correlations may indeed have a certain significance. 

Quantitative understanding of the measured excitation 
spectrum and of the underlying strong correlations in the 
2e QD were obtained through exact diagonalization of the 
Hamiltonian for two 2D interacting electrons 

U = H{v 1 ) + H{v 2 )+ 1 e 2 l{ K r 12 ), (1) 

where the last term is the Coulomb repulsion, k (12.5 for 
GaAs) is the dielectric constant, and r\i = |ri — v^\. The 
prefactor 7 accounts for the reduction of the Coulomb 
strength due to the finite thickness of the electron layer 
in the ^-direction and for any additional screening ef- 
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fects due to the gate electrons. H(r) is the single-particle 
Hamiltonian, which for an anisotropic, harmonically con- 
fined QD, is written as 

H(r)=T+^m*(u J 2 x x 2 +u J 2 y y 2 ), (2) 

where T = (p - eA/c) 2 /2m* , with A = 0.5(-By, Bx, 0) 
being the vector potential in the symmetric gauge. The 
effective mass is m* = 0.07mo, and p is the linear mo- 
mentum of the electron. The second term is the exter- 
nal confining potential. In the Hamiltonian J2J, we ne- 
glect the Zeeman contribution due to the negligible value 
(.9* ~ 0) of the effective Lande factor in our sample. 

In the EXD method, the many-body wave function is 
written as a linear superposition over the basis of non- 
interacting 2e determinants, i.e., 

2K 

^'L(ri,r 2 ) =J2 n « W 1 ; W(2;i)>, (3) 

where = (pi{v\)a{l) if 1 < % < K and tp(l;i) = 

<Pi-K(ri)/3(l) if K + 1 < i < 2K [and similarly for 
ip(2,j)], with a(P) denoting up (down) spins; the index i 
(or j) = (k,l) and <^(r) = X k (x)Yi(y), with X k (Yi) be- 
ing the eigenfunctions of the one-dimensional oscillator in 
the x(y) direction with frequency Lo x (u> y ). In the calcu- 
lations we use K = 79, yielding convergent results. The 
total energies -Eexd ano - the coefficients f2*j are obtained 
through a "brute force" diagonalization of the matrix 
eigenvalue equation corresponding to the Hamiltonian in 
Eq. (JTJ. 

Figure [3] shows the experimentally extracted energy 
splitting J(B) between the states So and T + together 
with results of the EXD calculations with huj x = 4.23 
meV, ftwy = 5.84 meV [i.e., fkuo = h^J (us 2 + iOy)/2 = 
5.1 meV], and 7 = 0.862 With this value of 7, 

the Wigner parameter Rw = 1-34. The anisotropy used 
in the calculations is Lu y /u> x = 1-38 indicating that the 
quantum dot considered here is closer to being circular 
than in other experimental systems 0, . 

To assess the dot parameters determined from the 
comparison between the measured and calculated J{B), 
we checked that varying the value of loq (while keeping 
ijJy/oJx — 1-38) and the parameter 7 does not influence 
substantially the quality of the fit to the experimental 
J(B) — as long as huj does not deviate greatly (i.e., 
more than 20%) from the N = 1 single-particle excita- 
tion Ai. Such range of variations in the dot parameters is 
reasonable for the given sample and its measured N = 1 
excitation and charging energies and is due to changes in 
the confinement with changing V pg . 

The inset of Fig. 01 shows the EXD calculated singlet 
conditional probability distributions (CPDs) [l6j, [TJ, ]2li 
I24I [25| for the second electron, given that the first elec- 
tron is kept fixed at a given location (indicated by a solid 
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FIG. 3: (Color online) The ST splitting J(B) extracted 
from the experiment (open squares) together with the cal- 
culated EXD results (solid line). For the dot parameters, 
see text. Inset: Conditional probability distributions for the 
second electron when the first is held fixed at the point in- 
dicated by the solid dot. For the 2e QD that we studied, 
the Wigner molecule character is gradually enhanced with in- 
creasing magnetic field. Lengths in nm. CPDs in arbitrary 
units. 

dot). The electrons clearly tend to avoid each other re- 
flecting their strongly correlated motion, in contrast to 
the two electrons in the natural Helium atom that oc- 
cupy the same single- particle orbital. Comparison of the 
CPDs for B = and B = 3.8 T shows that the mag- 
netic field leads to an enhancement of the correlation be- 
tween the two electrons resulting in a sharpening of the 
Wigner-molecular state. Larger values of B result in an 
effective dissociation of the WM in analogy to the 
dissociation of the natural H2 [3|, leading eventually to 
a vanishing of J(B), in agreement with the trend seen in 
Fig. We emphasize that the negligible Zeeman split- 
ting in our sample is crucial for the observation of this 
trend which should even lead to ST oscillations of the 
ground state spins at higher magnetic fields ■ Exper- 
imentally, however, the accessible magnetic field range is 
limited, because conductance resonances disappear with 
increasing field due to wave function shrinkage, and coun- 
terbalancing this effect by changing gate voltages would 
seriously change the confinement potential. 

The EXD calculated spectrum of the elliptic 2e QD is 
displayed in Fig. 4. The colored curves correspond to fea- 
tures highlighted in the experimental spectrum (see Fig. 
2) in the region between N = 1 and N = 2. The good 
overall agreement between the measured and calculated 
spectra allows for a full interpretation of the observed 
features and their physical origins rooted in strong elec- 
tronic correlations. The energy difference between the 
T + and Sq states and its variation with B were discussed 
above [see J(B) in Fig. 3]. 

Further interpretation is facilitated by consideration of 
the EXD calculated spectrum of the corresponding cir- 
cular dot (Fig. 4, inset) where all the many-body energy 
levels (0, m) with zero nodes and angular momentum in- 
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dex m > decrease in energy with increasing magnetic 
field, converging at much higher B to the lowest Landau 
level of the interacting electrons. Since the deformation 
of the elliptic QD discussed here is moderate, we adopt 
in Fig. 4 the same notation as for the circular dot (see 
caption of Fig. 4). The second excited singlet state seen 
in the experiment (Fig. 2) and its avoided crossing with 
the first excited singlet state at B ks 5.8 T is well repro- 
duced by the theory [see curves marked as S2 (red) and 
So (magenta) in Fig. 4; the avoided crossing is marked 
by a double arrow] |2g. The third excited state, T+ t cM, 
is identified as a (1,0,0,1) triplet excited state (green on- 
line) that is related to the T + state (blue online) through 
a CM excitation along the x direction. We note here that 
the aforementioned avoided crossing between the Sq and 
the 52 excited states developed as a result of the shape 
deformation of the dot and out of the crossing between 
the (0,0) and (0,2) singlets of the circular dot (marked by 
a circle in the inset). Another correspondence between 
the measured and calculated spectra is the splitting at 
B = between the T + and T_ states (measured as 1.5 
meV and calculated to be 1.9 meV); these triplet states 
are degenerate for the circular dot (see inset). 
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FIG. 4: (Color online) Calculated EXD energy spectrum, ref- 
erenced to 2h^/u}Q + oi§/4, in a magnetic field of a 2e dot 
with anisotropic harmonic confinement (for the dot parame- 
ters, see text). We have adopted the notation (N x , N y ,n,m), 
where (N x ,N y ) refer to the CM motion along the x and y 
axes and (n, m) refer to the number of radial nodes and an- 
gular momentum of the relative motion in the corresponding 
circular dot. Inset: The EXD spectrum of the correspond- 
ing circular dot. Only the (n, m) indices are shown, since 
N x = N y — for all the plotted curves. Solid lines denote 
singlets. Dashed lines denote triplets. 

In conclusion, we have presented measurements of the 
excitation spectra of a lateral 2e QD, with negligible 
Zeeman splitting, moderate anisotropy and parabolic z- 
confincmcnt being the main differences to systems in- 
vestigated before. Exact-diagonalization calculations, in- 
cluding a rather small (14%) decrease in the effective in- 
terelectron Coulomb repulsion, explain in detail the ob- 
served features of the excitation spectra under the influ- 



ence of a variable magnetic field, including the variation 
of the ST splitting, J, as a function of B. Theoretical 
analysis of the salient excitation-spectra patterns (e.g., 
avoided level crossings, the line shape of J(B), and the 
lifting of level degeneracies in comparison to the circular 
case) indicates significant electronic correlations, indicat- 
ing the formation of a well developed Wigner molecule in 
an anisotropic quantum dot. In particular, the behav- 
ior of J(B) for larger fields suggests an analogy between 
the dissociation process of the WM and that of the nat- 
ural H 2 molecule. These results are expected to assist 
the design and implementation of quantum-dot based de- 
vices, as well as to further advance our understanding of 
the interplay between the confinement strength, shape 
anisotropy, and degree of screening that also character- 
ize other recent measurements on few-electron quantum 
dots 
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